Confinement Induced Quantum Phase Transition and Polarization Cooling in a 

Dipolar Crystal of Polar Molecules 
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It is well-known that the liquid properties in a strongly confined system can be very different 
from their ordinary behaviors in an extended system, due to the competition between the thermal 
energy and the interaction energy. Here we show that, in a low- dimensional self-assembled dipolar 
crystal, the parabolic optical confinement potential can also strongly affect the quantum many-body 
properties in the low temperature regime. For example, by changing the confinement aspect ratio, 
the bulk of the system can undergo a quantum phase transition between a liquid state and a solid 
state via a nonmonotonic pattern formation of the domain wall. Furthermore, the entropy of a 
trapped dipolar crystal can be much larger than the liquid state in the weak dipole limit, indicating 
an intrinsic polarization cooling mechanism via increasing the external field. These highly correlated 
confinement effects are very important to the experimental preparation of a self-assembled dipolar 
crystal using ultracold polar molecules. 
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Confinement effects are known to be very crucial to 
the physical properties of a spatially confined liquid. In- 
vestigating these extraordinary fluid properties has be- 
come an important subject in nanoscience and technol- 
ogy in recent years [l[. When the confined length scale 
is reduced to be comparable to the inter-particle dis- 
tance, there could be several kinds of phase transitions 
observed in different systems, including solidification on 
the boundary Q, liquid-to-liquid phase transition Q, 
and layering phenomena [J], etc. From a microscopic 
point of view, all these unusual properties result from 
the competition between the thermal kinetic energy and 
the inter-particle interaction in a certain spatial geom- 
etry. It is therefore reasonable to expect that the in- 
terplay between quantum fluctuation effects, interaction 
energy, and the confinement geometry can also lead to 
some extraordinary many-body properties that are not 
observable in classical systems. Here we show that in a 
self-assembled low-dimensional dipolar crystal Q , which 
can be formed by ultracold polar molecules trapped in 
a magneto-optical potential, the bulk properties of the 
crystal can be strongly affected by the confinement po- 
tential, even when the system size is still much larger 
than the average inter-particle distance. The thermody- 
namical properties are also found qualitatively different 
from their behavior in an extensive space, showing the 
highly correlated confinement effect of a quantum many- 
body system. Different from the observed dipolar and 
confinement effects in the condensate of magnetic atoms 
Q , here we mainly focus on the dipolar crystal phase in 
the strongly interacting regime (see below), which can- 
not be easily achieved by magnetic atoms within present 
experimental parameter range. 

We first consider polar molecules initially prepared in 
the lowest rotational ground state (L = 0) and confined 
in a quasi-two-dimensional (2D) trap, where the trans- 
verse dynamics (in z axis) is frozen to the single particle 
ground state. The in-plane motion of molecules, however, 
is weakly confined by weaker harmonic potentials with 



trapping frequencies, oj x . yi in the x and y directions re- 
spectively. When an external DC electric field is applied 
along the z axis, polar molecules become polarized and 
have a field-dependent electric dipole moment, D. In a 
dilute limit, the mutual interaction between these quais- 
2D polar molecules is a dipolar interaction and the 
system can be easily stabilized by the strong transverse 
confinement potential. As a result, the system Hamilton- 
ain can be written to be 
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where m is the mass and N is the total number of 
molecules, and — are the in-plane momen- 

tum and position operators of the ith particles. There 
are three length scales in Eq. JT]): two are the oscillator 
lengths, <x OSC)X («) = y/W™^xM> an d one is the dipolar in- 
teraction strength, ad = mD 2 /h 2 . The system properties 
are therefore determined by three dimensionless param- 
eters only: 7 = ad/a osc ^ x , measuring the dipolar interac- 
tion strength, k = us y /uj x — aosc,i/ a osc,y > measuring the 
aspect ratio of the confinement potential, and particle 
number N. Throughout this paper, we use a osc = a osc . x 
and uj = uj x to be the units of length and energy scales 
respectively. For a typical molecule like SrO, the full po- 
larized dipole moment is D = 8.9 Debye and therefore 
7 ~ 88 for a osc — lAfim in a typical trapping frequency, 
uj = 2tt x 50 Hz. 

In this paper, we are interested in the regime when 
the dipolar interaction is sufficiently large (7 > 1), so 
that polar molecules can form a dipolar crystal [g], as 
shown in Fig. [Ha) (for ID trap, k = 00) and Fig. 
[TJe) (for 2D isotropic trap, k = 1). Here the equilib- 
rium position for each molecule is calculated by using 
Molecular-Dynamics (MD) simulation with a small fric- 
tion in the Langevin equation [8j|9j . We further calculate 
the phonon spectrum by quantizing the position fluctu- 
ation of dipoles from their equilibrium positions to the 
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quadratic order. Within this harmonic approximation, 
the obtained spectrum (Figs. QJb) and (f)) is indepen- 
dent of the dipolar strength (7) , because the effect of tun- 
ing dipole moment can be exactly cancelled by adjusting 
the inter-molecule distance [13 ■ In Fig. Qlc)-(d) and 
(g)-(h), we show some typical phonon excitation wave- 
functions for these two systems respectively. One easily 
see that the confinement potential has strongly changed 
the excitation wavefunctions from simple plane waves in 
a uniform system. 

Following the same spirit of the Lindermann criterion, 
which was first applied to the classical 3D melting prob- 
lem [12] and then extended to the 2D thermal and quan- 
tum melting [13lll4lll5l |. in this paper we define that a 
solid phase is melted at a certain position if its posi- 
tion fluctuation is larger than an imperical ratio (C'l) of 
the lattice constant. Such ratio in a uniform 2D dipolar 
crystal has been calculated to be 0.23 from the quantum 
Monte Carlo simulation (l^flB |. To apply such useful 
criterion to a nonuniform system discussed in this paper, 
here we define the local "softness" of the dipolar crystal 
to be: & ee An/k, where Ar l = yj ((r. ( - 170) 2 ) is the 
position fluctuation of the ith particle (calculated from 
the wavefunction of the phonon excitation states), and 
Zj is its average distance to the nearest neighboring sites. 
As a result, the ith particle is considered to be in a liquid 
state if > Cl , while it is in a solid state if & < Cl [3 • 
In Fig. [U(a) and (b), we show the local softness (£,•) for 
dipolar crystal in both ID and 2D isotropic potentials. 
We find that in ID system the softness distribution is 
convex with its maximum value in center of trap, i.e. the 
bulk of the ID dipolar crystal is softer than the edge, in 
contradiction to a naive guess from a local density ap- 
proximation, which predicts the position of higher parti- 
cle density should be more solid-like due to the r -3 dipo- 
lar interaction. This striking results originate from the 
confinement effect on the quantum fluctuation of lattice 
points: the system ground state wavefunction has equal 
contribution from each eigenmodes (note that quantum 
zero point energy is Eq = J2 n ^^n), while higher energy 
modes prefer to generate much larger position fluctuation 
in the center of the ID trap (see Fig. Q] (c) and (d)). In 
the 2D isotropic trap, on the other hand, the softeness be- 
comes concave, because now the high energy eigenmodes 
can still have large amplitude at the edge due to the ad- 
ditional degree of freedom in the azimuthal direction of 
the trap (see Fig. QJg) and (h)). In Fig. [2]Jc) and (d), we 
further show the phase boundary (domain wall) between 
the liquid state and the solid(crystal) state as a func- 
tion of dipolar strength (7): the ID system has a central 
domain of liquid state as long as 7 is smaller than a crit- 
ical value, 7i D , above which the whole system becomes 
a well-defined solid crystal. On the other hand, the cen- 
ter of the 2D isotropic trap starts to be crystalized only 
when 7 is larger than another critical value, Such 
qualitative difference in these two systems are clear evi- 
dence of the highly correlated confinement effects on the 
quantum fluctuations. Both 7* D and 72 D decrease as the 



number of particles increases. 

In Fig. [3l we show the confinement effects in different 
trap aspect ratio with a fixed lu x : the anisotropic trap- 
ping potential generates two local minimum of the soft- 
ness (see Fig. OJb)) along the x axis. Within a proper- 
parameter range, these two local minimum of softness can 
become the centers of two solid island(domain) embedded 
in a dipolar liquid. In Fig. [3^c), we show how the local 
softness, in the edge (red lines) and in the center (black 
lines) of the trap, changes as function of the aspect ratio, 
k. We find that although the former decreases monoton- 
ically as expected, the later can have several re-entrant 
effects and can be even larger (i.e. softer) than its value in 
the isotropic trap (n = 1) without compression. Increas- 
ing number of particles does not change these nonmono- 
tonic bulk properties (say domain walls and re-entrant 
behavior), showing a significant confinement effects even 
when the system size is much larger than inter-particle 
distance. Similar nonmonotonic/reentrant behavior are 
also observed in a system of classical 2D melting [lj| but 
much more particles and quantum fluctuation effects are 
considered here. 

The nontrivial confinement effect can be also observed 
in the finite (but low) temperature regime by investigating 
the system entropy (S) , which can be assumed to be con- 
served during an adiabatic manipulation of the system 
parameters. The total entropy can be calculated from 
0S = Jo dT'C v (T')/T', where C V (T) = dE(T)/dT 
is the specific heat. The total energy, E(T), can be 
easily calculted from the phonon excitation spectrum: 
E(T) =E C + E Q + J2 n js&frzi , where E c is classical 
potential energy and ks is Boltzman constant. In Fig. 
EKd), we show how the system temperature changes as a 
function of the confinement aspect ratio, k, by keeping 
the total entropy a constant during the adiabatic pro- 
cess: the temperature increases sublinearly as the sys- 
tem is compressed and eventually becomes saturated in 
the limit of pure ID system (say k > 100 for N = 91). 
Within the harmonic expansion approximation for the 
phonon excitations, the calculate system entropy is inde- 
pendent of the dipolar strength, 7. 

From experimental point of view, one of the most im- 
portant question is how the system temperature changes 
when a dipolar crystal is formed by increasing the elec- 
tric field. In the zero dipole moment limit (i.e. zero 
external field), only s-wave scattering exists between 
bosonic molecules, while it is almost noninteracting for 
single component fermionic molecules. For a Bose liq- 
uid in a 2D isotropic trap, we can apply the local den- 
sity approximation to calculate the system entropy [20| . 
and obtain the leading order temperature dependence: 
S B h/k B = a^L (^Z.) 2 with a ~ 5.7 x 10~ 3 being ob- 
tained from direct numerical calculation. This result is 
independent of the s-wave scattering length and parti- 
cle number due to the unique geometr y o f 2D isotropic 
harmonic trap. Similar to the 3D case [20j], it applies to 
condensate as well as to the normal state, because the 
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dominate contribution of entropy always comes from the 
single particle excitation of normal liquid. For the nonin- 
teracting fermionic molecules, we can also calculate the 
system entropy easily from the temperature dependence 
of the total energy in a 2D harmonic trap [2l| , and obtain 

S FL /k B = (£gf) . In Fig. |2 we show the calcu- 

lated entropies of Fermi liquid (dotted lines), Bose liquid 
(dash-dotted line) and dipolar lattice (solid lines) in a 2D 
isotropic trap. Results for both N — 91 and N = 217 
are shown together for comparison. It is easy to see that 
if bosonic polar molecules are initially prepared at zero 
field and in a sufficiently low temperature (< T£ OBon ), 
the system temperature will decrease (intrinsic cooling, 
blue leftward arrow) as the external field is increased to 
derive the system toward a dipolar crystal adiabatically. 
On the other hand, if the system is prepared in a rather 
high initial temperature (> Tb oson ), the system tempera- 
ture will increases greatly as the dipole moment increases 
(intrinsic heating, the red rightward arrow). The critical 
temperature, T£ oson ~ lQhuj/k B for N = 91, but be- 
comes about lQfvjj/ks for N = 217. This shows that 
the polarization cooling here is due to the many-body 
effects of polar molecules, completely different from the 
demagnetization cooling process in solid state systems or 
in magnetic dipolar atoms [22j. More precisely, due to 
the presence of a harmonic confinement, the average dis- 



tance between dipoles increases (i.e. the average density 
decreases) as the dipole moment is enhanced, and there- 
fore the system temperature can be reduced via transfer- 
ring the electric field energy to the confinement potential 
energy (rather than to the kinetic energy) in the limit of 
strong dipolar crystal. Such intrinsic cooling mechanism 
is totally different from what is expected in a uniform 
system, where the entropy of a crystal should be always 
smaller than a liquid state at the same density 23]. Sim- 
ilar intrinsic cooling mechanism can be also observed in 
fermionic polar molecules and/or in ID trapped system. 

In summary, we have shown several important confine- 
ment effects on a self- assembled dipolar lattice formed by 
ultracold polar molecules. Changing confinement aspect 
ratio can induce a quantum phase transition between a 
bulk liquid state in ID trap to a bulk solid crystal state 
in 2D isotropic trap. We further find an intrinsic polar- 
ization cooling mechanism during the formation of dipo- 
lar crystal. Our results can be applied to the experi- 
mental preparation of a dipolar crystal of ultracold polar 
molecules and are also important to the understanding of 
the confinement effects in a quantum many-body system. 
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B.I. Halperin, and Lin I. This work is supported by Na- 
tional Science Concil. 
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FIG. 1: (a) Density profile and position configuration of 
lattice position for N — 30 dipolar particles in a ID har- 
monic trap with 7 = 17. (b) Excitation spectrum for the 
same system with N = 91. The red line is obtained in 
the hydrodynamic theory of the liquid state with N — > 00: 
lu„ — uj x -^/(3n 2 — n)/2 (c) and (d) show the eigenfunc- 

tion (i.e. deviation from their equilibrium position) for the 
lowest (n = 1, 2, 3) and the highest (n = 91) excited states of 
ID system, (e)-(h) show the same physical quanities as (a)- 
(d), but for N = 91 particles in a 2D isotropic harmonic trap 
with 7 = 80. (g) and (h) are the second and the seventh exci- 
tation modes with vibrating direction shown by arrows. The 
red lines in (h) are eye-guiding, indicating a Tachenko's mode 
as observed in the vortex lattice of a fast rotating condensate 





FIG. 2: (a) and (b): The distribution of softeness for 
dipolar crystals with N — 91 particles in ID and 2D isotropic 
traps respectively. 7 = 80 in (b). (c) and (d): The radii of the 
phase boundary, R, as a function of dipolar strength, 7, for 
the above two systems with N = 91. iitotal is the radius of the 
whole system. Here we set Cl = 0.23, and phase boundary 
for JV = 217 are also shown for comparison. 
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FIG. 3: (a) Position configuration of the dipolar crystal 
(N — 91) in elliptic traps of different aspect ratios, k. (b) 
Distribution of the softness for k — 2. (c) Local softness in 
the center (filled squares) and at the end (open squares) for 
systems of different numbers of particles (N). (d) The system 
temperature changes during an adiabatic compression from a 
2D isotropic trap (k = 1) to ID trap (k 3> 1) with N = 91. 
The confinement frequency in the x direction is fixed, and all 
results of are calculated with 7 = 80. 




FIG. 4: Total entropy as a function of system temperature in 
a 2D isotropic confinement potential. Results for Bose liquid 
(dash-dotted lines), Fermi liquid (dotted lines) and crystal 
phase (solid lines) are shown together. Results for different 
numbers of particles are also shown for comparison. The blue 
leftward arrow indicates the intrinsic cooling process for a 
Bose liquid, while the read rightward arrow indicates the in- 
trinsic heating process. 



